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ABSTRACT 

We have discovered a partially eclipsing white dwarf, low-mass M dwarf binary (3.015114 hour orbital 
period), SDSS J143547. 87+373338.5, from May 2007 observations at the WIYN telescope. Here we present 
blue band photometry of three eclipses. Eclipse fitting gives main sequence solutions to the M dwarf companion 
of M s = O.15-O.35M and R s = 0.17-0.32fl Q . Analysis of the SDSS spectrum constrains the M dwarf further 
to be of type M4-M6 with Ms = 0.11- 0.20M Q . Once full radial velocity curves are measured, high precision 
determinations of the masses and radii of both components will be easily obtained without any knowledge of 
stellar structure or evolution. ZZ Ceti pulsations from the white dwarf were not found at our 4 mmag detection 
limit. 

Subject headings: binaries: eclipsing — stars: white dwarf — stars: individual: SDSS J143547. 87+373338. 5 



1. INTRODUCTION 

Eclipsing binary systems offer unique opportunities for pre- 
cise measurements of stellar properties. With full radial ve- 
locity and light curves (commonly from eclipsing double- 
lined spectroscopic binaries), the masses and radii of both 
components are measured to very high precision and accu- 
racy. 

We discovered that the u 1 = 17.7 and g' = 17.1 
mag binary SPSS J1 43547.87+373338.5 (hereafter J1435, 
Eisen stein et alJl2006h undergoes a partial eclipse of the WD 
by its low-mass companion with a period of 3 hrs and a tran- 
sit time of w 480 sec. The eclipsing M dwarf secondary 
is expected to contribute < 4% in the blue band BG-39 fil- 
ter. Measurements by Gansicke et a l. (private communica- 
tion, Rebassa-Mansergas et al. 2007) give an accepted range 
of WD masses of 0.35 -O.58M (see Table Q] and ®. Our 
analysis in $3] using model fits to th e partial eclip s e prod uce 
main sequence (MS) solutions using Baraffe etal] d!998l) for 
the M dwarf of M s « 0.15-0.35M o and R s « O.17-O.32# , 
making it a candidate for probing the stellar structure of the 
low-mass MS. The M dwarf is further constrained in i2]by us- 
ing the SDSS spectrum and is expected to be of type M4-M6 
with Ms « 0.1 1 — O.2OM . Measuring the ingress and egress 
as well as ellipsoidal and reflection effects in multiple color 
filters of the J 1435 eclipse would probe the WD atmosphere. 

In $2]we detail our observations and data reduction. A deep 
analysis for ZZ Ceti type pulsations revealed no pulsations 
down to a 4 mmag detection limit. In ^3]we discuss our anal- 
ysis of the eclipse light curve and use it to constrain the prop- 
erties of J1435. In |2]we discuss the consequences our results 
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TABLE 1 

Properties of SDSS J143547. 87+373338 



Property 


Eisenstein et 


Gansicke et 


Gansicke et 




al. (2006) 


al. (Hot) 


al. (Cold) 


7k (K) 


11062 ±58 


12681 ±990 


12392 ±350 


logg (dex) 


6.860 ±0.038 


7.68 ±0.20 


7.82±0.14 


M W D (Mq) 


0.20 a 


0.35 b 


0.58 c 


Rwd (fi Q ) 


0.0273 a 


0.0178 b 


0.0132 c 



a From lAlthaus & Benvenuto 11997) He core model , 
using logg = 6.86. b From lAlthaus & Benvenuto! 09971) 
He core model, using logg = 7.48 lower limit. From 
lAlthaus & Benvenuto] ©98) C/O core with IO^Mq H enve- 
lope model, using logg = 7.96 upper limit. 

along with the SDSS spectrum have on the properties of the 
secondary star. 

2. OBSERVATIONS AND DATA REDUCTION 

SDSS J 1435 was targeted, along with many other objects, 
in a campaign to discover very-low-mass (He core) and very- 
high-mass (O-Ne core) pulsating DA WDs (ZZ Cetis) . From 
within the Sloan Digital Sky Survey, Eisens tein et al .1 d2006l) 
released over 9,000 spectroscopically classified WDs. We se- 
lected only tho se DA objects within the empirical ZZ Ceti 
instability strip (Gianninas et al. 2006; Mukadam et al. 2004) 
of either low (log g < 7.45) or high (log g > 8.7) gravity. 
J 1435 was among several dozen objects that met the low grav- 
ity criteria and was visible during our observing run of 29 
May through 1 June 2007 . In Table [T| we l ist the properties 
of J1435 as measured by Eisenstein et al. (2006) and Gan- 
sicke et al. (priv. com., Rebassa-Mansergas et al. 2007). We 
report both a "cold" and "hot" solution for the WD parame- 
ters as the equivalent widths of the Balmer lines go through 
a maximum near r e ff « 13,000 K (dependent on logg) and 
two solutions of similar quality result on either side of this 
maximum. In most cases, this degeneracy is lifted by fitting 
the overall shape of the spectrum, however, since J1435's T e ff 
is so near the maximum, both solutions are indistinguishable 
(Rebassa-Mansergas et al. 200 71). The large discrepancy in 
parameters is because lEisenstein et all d2006l) did not spec- 
trally subtract the comp anion's contribution whereas Gan- 
sicke et al. (priv. com., Rebassa-Mansergas et al. 2007) did. 

We observed J 1435 on the nights of 29 May, 31 May, and 1 
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FIG. 1. — Top Panels: Light curves of all three nights in relative magni- 
tudes. Bottom Panel: Phased light curves of all nights centered on the pri- 
mary eclipse and in relative flux units. Note the representative average error 
bar of the measurements (grey points). The solid line is the eclipse fit (see 
i(3) for Mwd = 0.35Mq (first line in Tableff). The other two lines are the best 
possible fits for other MS masses at Mwd = O.35M0 . 



June 2007 with the 3.5-meter WIYN telescope at the Kitt Peak 
National Observatory. All observations used the OPTIC cam- 
era with two 4Kx2K pixel CCDs side-by-side (15 - /im pix- 
els) for a total of 4Kx4K pixel viewi ng area with a full frame 
field of view of 9.5 x 9.5 arcminutes (iHowell et alJl2003l) . Ex- 
posure times were 15 seconds over w 2 hr for the 29 May 
and 31 May observations, and 40 seconds over sa 3.5 hrs on 
1 June. All exposures were with the broadband BG-39 fil- 
ter (A c « 4800A, FWHM « 2600A). The 4Kx4K total pixel 
CCD was binned 2x2 to reduce the readout time to w 8.1 
seconds, which varied by » 0.03 seconds during a night. 

All images were reduced by applying an averaged two di- 
mensional bias subtraction, and were flattened using averaged 
and normalized dome flat field images taken immediately 
prior to the night's observations. All tasks were performed us- 
ing standard tasks within IRAF 6 . Light curves were produced 
by extracting fluxes from the program star and several com- 
parison stars and taking their rati os. Fluxes were extra cted us- 
ing the IRAF package VAPHOT dDeeg & Dovld2001h . which 
applies a dynamic method to aperture size adjustment (both 
flux apertures and sky background estimation annuli) to ac- 
count for variable seeing. Fifteen comparison stars were se- 
lected within the frame ranging from B ss 16 to 19.5. Differ- 



b c ompa 

weighting scheme detailed in ISokoloski et ail d2001l) and in 



ential photometry was performed usin g the c omparison star 

„■ n g I 

spired bv lGilhland & Brownl d!988l) . 

Though targeted as a possible ZZ Ceti using the r e ff 
and \ogg measured by lEisenstein et all (120061) . the new 7^ 
and log,? measurements by G ansicke et al. (priv. com., 
IRebassa-Mansergas et alj 120071) place it outside the empir- 

6 IRAF (Image Reduction and Analysis Facility) is distributed by the Na- 
tional Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under contract with the Na- 
tional Science Foundation, http://iraf.noao.edu 



ical i nstability strip dGianninas et alj 120061 : iMukadam et alj 
2004). It is therefore not surprising that our Lomb-Scargle 
periodogram analysis found no (to 4 mmags) pulsations in the 
frequency range of 1-17 mHz. 

3. ECLIPSE ANALYSIS 

All three nights of observation show distinct primary 
eclipses (see Figure [T|). Us i ng th e technique devel- 
oped by iKwee & van Woerdenl d!956h . fits to the tran- 
sit centers of these eclipses show primary minima at 
HJD 2,454,249.711028 ±0.000016, 2,454,251.7211338± 
0.0000044, and 2, 454,252.8517916± 0.0000051. Com- 
bining these primary minimum epochs gives an ephemeris 
of HJD 2,454, 249.71 1056 ±0.000011 + N 0.12562974± 

0. 00000055. Our analysis obtains an orbital period of P = 
3.0151 14±0.000013hrs (10, 854.410±0.048 sec) and a tran- 
sit time (measured as first departure from and full recovery to 
maximum light) of r w 480 sec. 

We constrain the secondary through detailed modeling of 
the eclipse light curve. There is no evidence of a secondary 
eclipse, so we only model the primary eclipse and assume 
the secondary star to be a black (non-luminous) disk. If we 
assume a point light source WD, Rwd *C Rs, Kepler's Law 
and the system's partially eclipsing geometry lead to a de- 
generacy between Ms, Rs, and i. We can connect Ms and Rs 
with a MS M-R relation, but are unable to break the rela- 
tionship to i. This last degeneracy is broken by giving the 
WD its finite physical extent, Rwd, and modeli ng the WD 
disk w ith a standard linear limb darkening law dVan H amme 
[19931) : 7(60/7(0°) = 1-m ld [1-cos(6I)], where is the angle 
made by the normal vector to the WD surface to the line- 
of-sight. We obtain the range M WD = 0.35-0.58M Q by us- 
ing the logg and r e ff ranges measured by Gansicke et al. 
(priv. com.. IRebassa-Mansergas et al.l 120071 see Table Q]). 
The low mass limit uses the He core WD cooling tables of 
lAlthaus & Benvenutol d!997[) while the high mass limit uses 
the cooling tables of lAlthaus & Benvenutol d!998l) for a C/O 
core WD with a IO^Mq H envelope of solar metallicity. Fi- 
nally, we assume circular orbits and the geometry of an in- 
clined system of orbiting masses. In all, there are seven pa- 
rameters: Mwd, Ms, Rwd, Rs, h P, and uio- We fit simple 
cos($) and cos(2$) functions to determine the strength of re- 
flection and ellipsoidal variations (respectively) on the light 
curve, but found no such variations at the 0.5% level. We 
neglect the modeling of these variations. 

In Figure [2] we present the three-sigma x| t contours for the 
Rs - Ms parameter space for two values of the linear limb 
darkening coefficient and extreme WD masses of the accepted 
range. We obtain a contour by fixing P, uio, Mwd, and Rwd- 
Then at each point in the Rs - Ms parameter space we min- 
imize %| t with respect to the inclination and adopt this x| t 
value for that mass and radius. Our data tightly constrains the 
radius of the secondary while very loosely constraining the 
mass if we assume a specific linear limb darkening coefficient 
and WD solution with non-negligible radius. Figure Q] shows 
a comparison of model fits to the observed light curve for one 
set of WD parameters and three sets of MS parameters fit over 

1. This figure plainly shows the broken degeneracy between 
M s -Rs and i. 

When modeling the limb darkening of a WD atmosphere 
with a linear law, the coefficient depends upon the effective 
temperature and the bandpass of the filter. Insufficient data 
exists to define this dependence, so we used two values brack- 
eting the measured range of similar systems: uid = 0.3-0.5 
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Fig. 2. — Constraints within the Rs — M$ parameter space. Three-sigma 
contours of minimized xL to eclipse light curve (see © using the two ex- 
treme accepted WD masses and t wo linear l i mb da rkening coefficients. RL 
overflow boundary defined using Eggleton ( 1983). Theoretical MS from 
Bara ffe et alj (\99S) . 
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FIG. 3. — Upper Panel: Constraints within the M$—Mwd param- 
eter space. All curves incorporate a minimized \ 2 fit to eclipse light 
curve (see j|3). All WD M-R r elations are derived as follows: He WD 
from[Al thaus & BenvenutcJ 11 9971) and C/O WD fr om Althau s"&~Benvenutd 
(1993) . The MS M -R relation is fro m lBaraffe et aT] 419981) . The RL bound- 
ary is defined using Eggleton 1 1983). Region I requires a secondary which 
is denser than a MS star of its mass, a likely sign of He enrichment. Region 
II requires a secondary which is less dense than a MS star of its mass, most 
likely not in thermal equilibrium and still contracting to the MS. Lower Panel: 
Radial velocity amplitudes for the WD and secondary for the MS solutions 
in the Upper Panel. 



TABLE 2 

Best Fit Main Sequence Secondary Solutions 



M W d 


ULD 


M s 


Rs 


i 


Vr,wd 


Vr, s 


(M ) 




(M Q ) 


(Re) 


(deg) 


(km/sec) 


(km/sec) 


0.35 


0.3 


0.143 


0.167 


79.3 


95.7 


234 


0.35 


0.5 


0.148 


0.172 


79.0 


98.3 


233 


0.58 


0.3 


0.340 


0.320 


72.4 


146 


248 


0.58 


0.5 


0.349 


0.327 


72.0 


148 


246 



dLittlefair et al.ll2006l l2007t iMaxted et alJl2007h . 

Given the present uncertainties and range of RL fillings, 
we simply assume that th e secon dary star has the Ms - Rs re- 
lation from iBaraffe et all (119981) . The MS solutions for the 
two extreme accepted WD masses and two limb darkening 
coefficients are in Table [2] Further, we investigated how the 
secondary star solution would change with the WD mass. In 
Figure [3] we plot the best fit M S star solutions using He core 
and C /O core WD models from lAlthaus & Benvenutol(ll997l 
1998) in the Ms — Mwd parameter space. 

We bound this space for low Ms by requiring that our sys- 
tem not be in RL overflow. If we assume our system to be at 
the boundary of RL overflow we can minimize xi t a l° n g the 
RL overflow boundary shown in Figure [2] The result of this 
minimization can be seen in Figure [3] as the boundary of the 
shaded region. The scaling of this curve is simply understood 
by considering the geometry of inclination and the circular 
orbits of two masses. By requiring an inclination where the 
WD is just fully eclipsed at zero phase, using the period and 
transit time of the binary, the radius of the secondary is found 
as a function of M wp and Ms. Equ ating this radius with the 
RL filling radius of Eggleton ( 1983) gives Ms as a function of 
Mwd whose scaling closely follows that of the curve in Fig- 
ure [3] Any system below this line must be in RL overflow. 
Figure[3]clearly illustrates two regimes for the secondary star. 
Region I requires it to be more dense than that of a MS star of 
the same mass, therefore, if found in this region it would be 
a likely sign of He enr ichment in the secondary f rom a prior 
stage of mass transfer dPylyser & Savonije1 fT988). Region II 
requires the secondary star to be less dense than that of a MS 
star of the same mass, therefore, if found in this region, the 
secondary would likely not be in thermal equilibrium. Both of 
these conditions are only possible if J 1435 is a mass transfer 
system that is temporarily out of contact. That being said, the 
most likely explanation for J 1435 is a pre-Cataclysmic Vari- 
able (pre-CV) where the M dwarf is a MS star that has yet to 
make contact. 

4. DISCUSSION 

Our analysis of J 1435 is far from complete. With a single 
color light curve of the primary eclipse, only broad constraints 
may be placed upon the parameters of the system (see ®. 
Therefore, we look to future observations in multiple colors as 
well as full spectra over the entire orbit to increase the preci- 
sion and accuracy of the system parameters. Our observations 
in a very blue bandpass, while ideal for a ZZ Ceti pulsation 
search, puts the effects of an M dwarf secondary eclipse at a 
level far below the precision of our observations. A redder 
bandpass would allow for better observation of the secondary 
eclipse and resulting model fits would provide additional con- 
straint to J1435's components. Full spectral coverage over 
the entire orbital period would ideally reveal a double-lined 
binary system yielding velocity curves for both components. 
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With these velocity curves, good measurement of the inclina- 
tion of the orbit can be made. 

One spectrum of J1435 was observed in the SDSS 
dYork et al.ll2000l) and within the errors of our ephemeris the 
entire exposure was taken out of eclipse. Visual inspection 
of the S PSS spectrum and co mparison with M dwarf spec- 
tra from iBochanski et al.l d2007l) clearly place the secondary 
star's spectral type to be later than M4, consist ent with under - 
filling the RL at P = 3 hrs with M s < 0.2M m dKniggef2 006). 
Using the empirical M dwarf colors from Boc hanski et al.l 
(2007) we fix the colors of the secondary at u' — g' > 2.3 and 
g' — i' > 3.0. Therefore, the M dwarf contribution in the BG- 
39 filter would be < 4% of the total system flux. This is con- 
sistent with the 10% remaining flux observed at maximum 
eclipse and the lack of an observed secondary eclipse. 

Further constraint of the secondary mass is possible by 
using the SDSS spectrum, synthetic WD spectra, and WD 
and M dwarf absolute magnitudes. We used a synthetic 
WD spectrum provided by G ansicke et al. (priv. com., 
Rebassa-Mansergas et alJ 120071) within the J1435 T e s and 
logg parameters and scaled it to fit the Balmer absorption 
lines shorter than 500 nm in the SDSS spectrum. In the John- 
son V filter at 560 nm, the contribution from a WD with our 
synthetic spectrum requi res the M dwarf to be » 2.3 mag- 
nitu des faint e r. Us ing iHolberg & Bergeronl d2006l) for WD 
and [Knigge (2006) for M dwarf absolute magnitudes, we 



find the M dwarf to be in the spectral range M4-M6 with 
Ms = 0.11- . 20M(T) . This analysis is in contrast to that of 
Silves trTet al.l {2006) which found the M dwarf spectral type 
of J1435 to be M2-M4 with M s = 0.20-0.55M Q . It is clear 
that much more rigorous spectral analysis needs to be carried 
out if it is to be used to constrain the secondary mass. 

We expect that J1435 is a pre-CV system with the M 
dwarf secondary on or near the MS. Two systems with sim- 
ilar parameters to J 1435, MS Peg an d NN Ser, have calcu- 
lated times to contact of w 1.5 Gyrs dSchreiber & Gansickd 
2003). There are little more than a dozen eclipsing double- 
lined spectroscopic binaries containing M dwarfs below 
\Mp, (j Bayless & Oroszll2 006: Heb b et alJl2006t iRibasi 120061: 
Young et al. 20 06tlDevor et al|2007 b lLo pez-Morales & Shaw 
l2007l: IShaw & L6pez-Moralesl2007l) . If J1435 were found to 
be a double-lined spectroscopic binary, from which velocity 
curves for both components could be found and its parame- 
ters found to high precision, it would be a significant addi- 
tion to this group especially being below 0.35M W , where the- 
ory dSiess. Forestini. & Dougados|[l997t iBaraffe et al.lll998t 
lYi et al. 2001) would benefit. 



We thank both referees for comments that clarified our pre- 
sentation. This work was supported by the National Science 
Foundation under grants PHY 05-51 164 and AST 07-07633. 
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